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Abstract 

We report on a search for ultra-high-energy (UHE) neutrinos from gamma-ray bursts (GRBs) in the data set col¬ 
lected by the Testbed station of the Askaryan Radio Array (ARA) in 2011 and 2012. Erom 57 selected GRBs, we 
observed no events that survive our cuts, which is consistent with 0.12 expected background events. Using NeuCosmA 
as a numerical GRB reference emission model, we estimate upper limits on the prompt UHE GRB neutrino fluence 
and quasi-diffuse flux from 10^ to 10*** GeV. This is the first limit on the prompt UHE GRB neutrino quasi-diffuse 
flux above 10^ GeV. 
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1. Introduction 

Gamma-ray bursts (GRBs) are the most powerful ex¬ 
plosions in the Universe. They emit high-energy gamma 
rays that are observable on Earth up to energies of ~ 100 
GeV, and are candidate sources of ultra-high-energy 
cosmic rays (UHECRs, above ~ 10^ GeV), whose ori¬ 
gin remains a mystery, and of neutrinos. The detection 
of neutrinos from GRBs would shine light on the ability 
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of GRBs to accelerate hadrons to the highest energies, 
and therefore on the possibility that they are the sources 
of the observed UHECRs. 

The widely accepted phenomenological interpreta¬ 
tion of particle acceleration in GRBs is the fireball 
model III] m a a 13. In this model, the energy car¬ 
ried by the electrons and hadrons in a jet of relativis¬ 
tic, expanding plasma wind — the fireball — may be 
dissipated through internal shocks between regions of 
plasma overdensity Elia. These shocks convert a sub¬ 
stantial part of the kinetic energy to internal energy by 
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accelerating the particles in the plasma. Accelerated 
electrons dissipate the internal energy as prompt gamma 
rays from synchrotron and inverse-Compton emission. 
Accelerated protons may dissipate the internal energy 
by interacting with the prompt gamma rays and produc¬ 
ing neutrinos in the 10^-10'° GeV range via a number 
of intermediate resonances nail. Later — typically, 
a few minutes after the prompt phase — the fireball 
collides with its surrounding medium, giving rise to re¬ 
verse and forward shocks. The latter are believed to be 
responsible for the GRB afterglow emission IfTOl [TTII . 
which may include neutrinos of energies comparable to 
the prompt ones ifT^ . 

Thus, GRBs might conceivably produce high-energy 
neutrinos copiously. However, due to the immense dis¬ 
tances separating us from the bursts — tens of Mpc 
to a few Gpc — the flux of neutrinos that arrives at 
Earth is expected to be low. Moreover, the flux is ex¬ 
pected to decrease with rising neutrino energy, due to 
the rising scarcity of protons of progressively higher 
energies at the sources. Over the last half-century, 
neutrino astronomy has steadily progressed in its abil¬ 
ity to detect low fluxes, culminating in the recent de¬ 
tection of a diffuse astrophysical neutrino flux, up to 
a few PeV, by the km-scale IceCube neutrino tele¬ 
scope unimiBKiiiniiisi. IceCube detects the op¬ 
tical Cherenkov light induced by neutrino interactions 
using > 5000 photomultipliers buried >1.5 km deep in 
the Antarctic ice. 

Significant sensitivity to higher neutrino energies re¬ 
quires larger detectors. While it can be cost-prohibitive 
to scale detectors that use techniques established for 
smaller scales up to volumes of order ~ 100 km^, an al¬ 
ternative is to utilize techniques that target a larger vol¬ 
ume with less instrumentation. 

One of the most promising methods to detect neu¬ 
trinos in the UHE range of 10*-10'° GeV in a large 
volume is the radio-Cherenkov technique Ifl9l . The in¬ 
teraction of a UHE neutrino in dense media induces 
an electromagnetic shower which develops a charge 
asymmetry. Because of this charge asymmetry, when 
the wavelength of the Cherenkov radiation is larger 
than the transverse size of the shower, the emission 
is coherent. This is known as the Askaryan effect 
Eol ED EH EH EH ESI Eg). For showers in ice, this 
process produces a radio frequency (RE) impulse at 
< 1 GHz which can be observed by antenna arrays 
read out with ~ GHz sampling rates. In this frequency 
range, the attenuation length in Antarctic ice is ~ 1 km 
Ell ESI, allowing a sparsely distributed array of de¬ 
tector units to observe volumes of ~ 100 km^. This 
is the strategy adopted by the Askaryan Radio Array 


(ARA) EDESEQ). In contrast, detectors that use op¬ 
tical Cherenkov signals are restricted by the < 100 m 
lengths over which attenuation, absorption, and scatter¬ 
ing diminish the signal, and thus require many more de¬ 
tector units to instrument the same volume ED. 

In this paper, we report on a search for UHE neutri¬ 
nos from GRBs using the 2011-2012 data set collected 
by the ARA Testbed station. Previous experiments have 
searched for neutrinos from GRBs using different tech¬ 
niques. However, they have either been sensitive to 
lower energies ll^ ES or only reported limits on the 
individual fluences of a handful of bursts EH- Instead, 
we present an upper limit on the stacked fluence of UHE 
prompt neutrinos from 57 selected GRBs and the first 
limit on the prompt UHE GRB quasi-diffuse neutrino 
flux in the range 10^-10^° GeV. 

This paper is organized as follows. In Section]^ we 
summarize previous GRB neutrino searches. In Sec¬ 
tion we describe ARA and the Testbed station. In 
Section |H we introduce our reference GRB emission 
model, NeuCosmA, and the AraSim detector simula¬ 
tion. In Section]^ we detail our data analysis pipeline. 
In Section]^ we present our results. In Section [7] we 
postulate future detection and analysis improvements. 
We conclude in SectionfS] 

2. Previous GRB Neutrino Analyses 

There have been many complementary GRB neu¬ 
trino searches reported by IceCube Il32ll35l[36l[37ll38l . 
ANTARES E3|39], RICE gO), and ANITA ll3Hl . 

IceCube ED is an in-ice, ~ 1 km^ optical-Cherenkov 
detector located at the South Pole. It has reported the 
most stringent limit on the GRB quasi-diffuse neutrino 
flux from 10^ to 10^ GeV ||36l . IceCube initially used 
an analytical GRB neutrino model by Guetta et al. El, 
based on the Waxman-Bahcall (WB) model EH, but 
now uses a numerical flux calculation EHEa that is 
compatible with the one used in the present analysis, 
NeuCosmA B3l . 

ANTARES iMl is an optical-Cherenkov detector, 
similar to IceCube, but located in the Mediterranean 
Sea, and instrumenting a volume of only ~ 0.03 km^. It 
is sensitive to a similar range of neutrino energies as Ice- 
Cube. The latest GRB neutrino analysis by ANTARES 
was based on NeuCosmA; its GRB neutrino flux limit 
is approximately an order of magnitude weaker than the 
limit from IceCube E9l . 

RICE ll45l was an in-ice radio-Cherenkov detector 
located in the South Pole, operational until 2011, that 
instrumented a volume of ~ 25 km^. The GRB neu- 
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ARA Testbed Station 
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Figure 1: Schematic of the ARA Testbed station. The 
borehole numbers are indicated next to their locations. 
Boreholes 1 through 3 and Borehole 5 each have a pair 
of Vpol and Hpol antennas while Borehole 6 has two 
Hpol antennas (Borehole 4 was not filled). The maxi¬ 
mum depth of the borehole antennas is ~30 m. 

trino analysis by RICE was based on an analytical neu¬ 
trino flux model and set individual fluence limits on five 
GRBs, from 5 x 10^ to 5 x 10*^ GeV BOll . 

ANITA 1461 is a balloon-borne Antarctic experi¬ 
ment that has flown three times under the NASA long- 
duration balloon program, searching for neutrinos us¬ 
ing the radio-Cherenkov technique. From an altitude 
of ~ 37 km, ANITA can monitor an extremely large 
volume of Antarctic ice, ~ 1.6 x 10^ km^ 1471 . The 
ANITA GRB neutrino analysis l3^ was based on the 
analytic WB GRB neutrino flux model l42l and set flu¬ 
ence limits for 12 individual GRBs that occurred in low- 
background analyzable time periods during its 31-day 
flight. ANITA provided the most recent GRB neutrino 
fluence limit from 10*^ to 10'^ GeV. The limited livetime 
of a balloon experiment constrains the maximum num¬ 
ber of analyzable GRBs for ANITA and thus they could 
not set a quasi-diffuse flux limit, but instead set fluence 
limits for each individual GRB. 

3. The ARA Instrument 

The full proposed ARA detector, ARA37, would con¬ 
sist of 37 stations spaced 2 km apart at a depth of 
200 m. The first three design ARA stations (Al, A2, 
A3) were deployed in the 2011-2012 and 2012-2013 
seasons, while a prototype Testbed station, which we 
used for this GRB neutrino search, was deployed in the 
2010-2011 season. 


Figure[2shows the layout of the Testbed with the po¬ 
sitions of the five boreholes. Boreholes 1 through 3 and 
Borehole 5 each contain a pair of antennas consisting 
of one vertically polarized (Vpol) bicone antenna and 
one horizontally polarized (Hpol) bowtie-slotted cylin¬ 
der antenna. Borehole 6, instead, has two Hpol quad- 
slotted cylinder (QSC) antennas which were deployed in 
the Testbed to test the antenna design before deploying 
them in the deep stations. All borehole antennas have 
bandwidths from 150 MHz to 850 GHz. For the trigger 
and data analysis, we utilized only antennas in Bore¬ 
holes 1-3 and 5. The maximum depth of the borehole 
antennas in the Testbed is approximately 30 m. There 
are also three calibration pulser VPol and HPol antenna 
pairs that were installed at a distance of ~30 m from the 
center of the Testbed array to provide in situ timing cali¬ 
bration and other valuable cross checks related to simu¬ 
lations and analysis. A more detailed description of the 
Testbed station is in Refs. W\W1- 

4. Analysis Tools 

In order to estimate the expected GRB neutrino spec¬ 
tra, we use the NeuCosmA GRB neutrino model. In 
order to estimate the efficiency of the ARA Testbed, 
we use AraSim, the ARA detector simulation software. 
Highlights of NeuCosmA and AraSim are described in 
the following sections. 

4.1. GRB Neutrino Model: NeuCosmA 

NeuCosmA Il48l |43]| is a state-of-the art computer 
code to calculate the neutrino fluence from cosmic ac¬ 
celerators such as GRBs. It performs detailed and fast 
computation of neutrino production in photohadronic 
py interactions, via A-resonance, higher resonances, 
decay channels, multi-pion processes, and direct pro¬ 
duction modes, and includes energy-loss processes of 
the secondaries and neutrino flavor oscillations during 
propagation to Earth. NeuCosmA provides fast calcu¬ 
lation of neutrino yields beyond simple analytical esti¬ 
mates, which are typically limited in the number of pro¬ 
duction modes. For each GRB, it provides the energy- 
dependent flavor composition of the neutrino fluence at 
Earth, i.e., the ratio of each flavor to the total fluence, 
ife,® • fn,® • /r.e) ■ 

We use NeuCosmA with model parameter values in¬ 
ferred from the observed gamma-ray signal of a GRB 
to calculate its neutrino spectrum. These parameters 
are Tgo (the time in which 90% of the gamma-ray flu¬ 
ence is collected), a and fi (spectral indices of the Band 
function 1491 at low and high energies), Fpeak (the peak 
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energy of the gamma-ray spectrum), Fy (the integrated 
gamma-ray fluence), Emin and fimax (the minimum and 
maximum energy of the fluence), and z (redshift). We 
extract parameter values from the Gamma-ray Coordi¬ 
nates Network (GCN) catalog 15011511. For unmea¬ 
sured parameters, we use their default values from the 
GRB-web database 1521 1^ . For all GRBs, we assume 
that the bulk Lorentz factor of the fireball F = 316, 
the energy in electrons and photons is equal to the en¬ 
ergy in magnetic fields, and the ratio of energy in pro¬ 
tons to energy in electrons — the baryonic loading — 
fp - 10 l32]|43. These are the same choices as in pre¬ 
vious analyses 1^ [Ml lJ7l l39l 1^ . 

Synchrotron energy losses of secondary tt^, tt^, n^, 
and /I* in the magnetic field of the source 153] |5l af¬ 
fect the shape and flavor composition of the neutrino 
fluence l55l . The onset of synchrotron losses for muons, 
pions, and kaons, at progressively higher energies, leads 
to GRB neutrino spectra that, in general, exhibit three 
distinctive kinks; see curves for individual bursts in Fig. 
1^ These effects, together with the energy dependence of 
the proton mean free path and the interaction of protons 
with the full photon spectrum, result in a quasi-diffuse 
neutrino flux — the “numerical fireball calculation” in 
Ref. EH — that is up to one order of magnitude smaller 
than the analytical estimates 19) used in the first IceCube 
GRB neutrino search 1321 . 

Contributions from different modes are performed via 
“response functions,” which contain the relevant kine¬ 
matics, multiplicities, and cross sections, encoded in 
fast-access look-up tables. This method is fast and ac¬ 
curate up to PeV energies. At higher energies, rele¬ 
vant for the present analysis, this approach has prob¬ 
lems treating the rising complexity in interaction final 
states, and QCD-based Monte Carlo methods like those 
implemented in SOPHIA ll56ll would give more accu¬ 
rate results. However, we expect that the impact of the 
particle-physics uncertainties is smaller than that com¬ 
ing from ambiguities in the astrophysical modeling of 
GRBs, even after reduction of errors due to averaging 
over the distribution of astrophysical parameter values. 
We discuss these effects more below. We use Neu- 
CostnA in the entire energy range of our analysis to 
obtain limits that are methodologically comparable to 
those found by other experiments. 

Our neutrino production model assumes that pro¬ 
tons are perfectly confined by the magnetic held at the 
source, and that only the neutrons produced in py inter¬ 
actions contribute to the flux of UHECRs. This “neu¬ 
tron model” results in a strong correspondence between 
the UHECR flux and the neutrino flux, which is in ten¬ 
sion with the non-observation of neutrinos from GRBs 


by IceCube lIMlIMlIMllJTl . All previous GRB neutrino 
searches have assumed the neutron model, so we adopt 
it to allow direct comparison of our results to theirs. We 
have not considered neutrino production models where 
protons can leak out of the source without interacting. 
They can yield neutrino fluxes lower by as much as an 
order of magnitude liSTlfSSll . So can models where mul¬ 
tiple shell collisions occur in the jet, each one with dif¬ 
ferent emission parameters li59ll60ll6TI . 

4.2. Detector simulation: AraSim 

AraSim li^ is a Monte-Carlo simulation software 
package used within the ARA Collaboration to simu¬ 
late neutrino signals as they would be observed by the 
detector. It simulates the full chain of neutrino events, 
such as the passage of the neutrino through the Earth, 
radio-Cherenkov emission, the path and response of the 
emitted signal in the ice, and the trigger and data acqui¬ 
sition mechanisms of the detector, as described below. 

AraSim was used in this search to model the neutrino 
interactions and detector response in the same manner 
that it was used in the ARA Testbed diffuse search, 
but we provide relevant details here for completeness. 
AraSim generates neutrino events with uniformly dis¬ 
tributed neutrino directions and interaction point loca¬ 
tions chosen with a uniform density in the ice. At each 
energy, we take the average flavor ratio of all GRBs 
given by NeuCosmA, weighted by their relative fluence. 
To properly account for the directional dependence of 
the sensitivity, the event is weighted by the probabil¬ 
ity that the neutrino survived its passage through the 
Earth and reached the interaction point. Once a neutrino 
interaction location is chosen in the ice, an in-ice ray 
tracing algorithm (RaySolver) derives multiple source- 
to-target ray-trace solutions giving signal arrival times. 
Erom each ray-trace solution, the radio-Cherenkov sig¬ 
nal, including a phase response, is then calculated with a 
custom parameterized radio-Cherenkov emission model 
inspired by Ref ll62l . The modeled signal is generated 
for both the hadronic and electromagnetic portions of 
the shower separately, as they have different character¬ 
istic shower profiles. We do not currently model the 
Landau-Pomeranchuk-Migdal (LPM) li6^ IMl l65l ef¬ 
fect in our RE emission model. Instead, we apply a cor¬ 
rection factor to the effective volume for each energy bin 
based on the impact of the LPM effect on the sensitivity, 
using the simpler RE emission model from Ref. lIMl . 

We then apply detector properties to the signal, such 
as antenna responses, amplifier and Alter responses, 
noise figure, and trigger mechanism. The antenna, am¬ 
plifier, and Alter responses are modeled based on simu¬ 
lation and measurements, while the noise figure and the 
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trigger mechanism are calibrated to the Testbed data. 
When a simulated event passes the trigger, the wave¬ 
forms are written into the same format as the data so 
that the simulated events can be analyzed with identical 
software. 

5. Data Analysis 

For this GRB neutrino search, we selected for anal¬ 
ysis only those GRBs that occurred during clean data- 
taking periods and in a region of the sky that is observ¬ 
able by our detector. After the GRBs are selected, we 
use the same selection criteria for the RF neutrino can¬ 
didate events as in the ARA diffuse neutrino search ll29l , 
but we search in a narrow time window around each 
GRB event, and thus we can loosen some cuts. We use a 
blinding technique that draws on both the ones used for 
the ARA diffuse neutrino search and the ANITA GRB 
neutrino analysis 041 . 

Our analysis consists of three stages. First, we use a 
10% subset from the full ARA Testbed data set for the 
preliminary background analysis. To estimate the back¬ 
ground, we use two 55-minute time windows on either 
side of each GRB event that excludes a 10-minute sig¬ 
nal window centered on that event. We optimize the 
cuts in the background analysis windows for the best 
expected limit in the signal windows. Second, we look 
at the number of events in the background analysis win¬ 
dows in the remaining 90% of the data set to check the 
consistency with the estimate based on the 10% subset. 
Third, we search for neutrino events in the signal win¬ 
dows in the entire (10%-i-90%) data set (note that the 
signal windows in the 10% set were not used for back¬ 
ground studies). 

5.7. GRB Selection 

We started with the 589 GRBs that occurred from 
January 2011 to December 2012 over the entire sky. For 
this analysis, we selected those that occurred during pe¬ 
riods of clean data-taking and that fell within the field 
of view of our detector. We used the IceCube GRB cat¬ 
alog im, which is based on the GCN lEollSI], to find 
GRBs during the time period of interest. 

From the 589 GRBs, we first rejected GRBs that 
failed the Effective Live time Cuts. The Effective Live¬ 
time Cuts consist of three cuts which require a low back¬ 
ground level and stable data-taking. The first cut is a 
simple time window cut which rejected GRBs that oc¬ 
curred during periods of high levels of activity at the 
South Pole station in the 2011 to 2013 seasons, in order 
to avoid strong anthropogenic backgrounds; for each 
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Eigure 2: Expected event spectrum from a simulated 
neutrino sample generated from the fluences of the 
257 GRBs that survived the Effective Livetime Cuts. 
Here we have applied the same analysis cuts that are 
used for the ARA diffuse neutrino search Il29ll . The 
ARA Testbed is most sensitive at ~ 10^ ^ GeV for these 
NeuCosmA-generated GRB neutrino fluences. 

year, we rejected GRBs that occurred from October 22'“^ 
to Eebruary 16*. The second cut requires that the data 
is not contaminated by any strong continuous wave¬ 
form (CW) source by rejecting any GRBs that occurred 
within an hour of any run where 10% or more events 
are highly correlated with each other. The third and fi¬ 
nal timing cut is a livetime cut which requires the detec¬ 
tor to be running and stably storing data within an hour 
of each GRB. The livetime represents the fraction of a 
second that the trigger was available. If there was any 
second when the livetime of the detector was lower than 
10% during the hour before or after a GRB, we reject 
that GRB from our analysis. After applying the Effec¬ 
tive Livetime Cuts, 257 GRBs survived from 224 days 
of analyzable period of data taking. 

To these surviving GRBs, we applied an additional 
cut which requires that the GRB should be included in 
the field of view of the Testbed. In order to define a 
field of view for the Testbed, we first found the energy 
bin which is the most sensitive to neutrinos from GRBs. 

Eigure l^is the expected event spectrum from the 257 
GRBs after applying analysis cuts that are used for the 
diffuse neutrino search Il29l . It shows that the Testbed 


5 




Figure 3: Effective volume of the ARA Testbed as a 
function of the zenith angle (6^) of the neutrino travel 
direction with a neutrino energy of 10^'^ GeV. The field 
of view is defined as the Full Width Half Maximum 
(FWHM) of the effective volume, which is -0.4 < 
cos0y < 0.05. This field of view covers ~ 20% of the 
sky. A vertically up-going neutrino has cos 6^-1. The 
shape of this distribution is described in the text. 


is most sensitive to NeuCosmA-generated neutrino flu- 
ences from these GRBs at ~ 10^'^ GeV. We used a sim¬ 
ulation set with the full range of incident angles of neu¬ 
trinos at 10^ ^ GeV, and obtained the effective volume as 
a function of neutrino direction. 

The effective volume 14^ is obtained for each energy 
bin and each neutrino direction bin by 


Feff = 


K, 


^thrown 


• lllggCICl 


( 1 ) 



Figure 4: The distribution map of 57 selected GRBs in 
Testbed local coordinates. The blue band in the map is 
the field-of-view cut range defined in Fig. Note that 
cos 6 in this map is the direction of the GRB while cos 6y 
in Fig.j^is the travel direction of the neutrino. 


Fig. which is -0.4 < cos Oy < 0.05. Earth absorption 
reduces the effective volume at high cos Oy (right-hand 
side of the plot), while the shadowing effect from the 
ray-tracing in ice causes the cut-off at low cos 0y (left- 
hand side of the plot) Il29l . 

Figure 1^ shows the distribution of the 57 GRBs that 
remain after applying a cut requiring that each GRB is 
within the field of view. They are shown in Testbed local 
coordinates, where </> - 0 points along the direction of 
ice flow and cos 0 = 0 points along the tangent to the 
surface of the geoid shape of the Earth. 

Figure 1^ shows the fluences of all 57 selected GRBs 
generated with the NeuCosmA software. Among them, 
one was brighter than the others; GRBl 10426A. Its flu- 
ence was higher than the others by an order of magni¬ 
tude for energies above 10’ GeV. Its location on the sky 
is marked as a cross in Fig.j^and its parameters values 
are shown in Table [T] The long duration and high spec¬ 
tral indices of GRB 110426A made its expected neutrino 
fluence significantly higher than for other GRBs at en¬ 
ergies above 10’ GeV. 

5.2. Neutrino search optimization 


where Vgen is a volume of ice where ice-neutrino inter¬ 
actions are generated uniformly, Vthrown is the total num¬ 
ber of events thrown (~ 10^ for each simulation set), and 
jnggered Weighted sum of the number of events 

that triggered. The weight Wi is the probability that the 
neutrino was not absorbed in the Earth, given its di¬ 
rection and the position of the interaction 

Figurej^shows the effective volume versus zenith an¬ 
gle of the neutrino travel direction. The field of view 
of the Testbed is defined as the Full Width Half Maxi¬ 
mum (FWHM) of the effective volume (arrow shown in 


This analysis uses the same set of cuts as in the Inter¬ 
ferometric Map Analysis in the ARA diffuse neutrino 
search ll29l . The analysis uses relative timing informa¬ 
tion to reconstruct the location of the source of the RF 
emission. The interferometric map is constructed from 
the sum of cross-correlations between the different pairs 
of antennas — a strong peak on the map indicates a high 
correlation among waveforms after correcting for the ar¬ 
rival times of the signals. We perform an optimization 
of the cuts for this analysis, which differs from the dif¬ 
fuse search by using the summed GRB fluence over the 
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GRB 

Tgo [sec] 

a 

P 

Fpeak [keV] 

Fy [erg cm 

Fmin [MeV] 

Fmax [MeV] 

z 

GRB 110426A 

376.05 

2.28 

3.28 

200 

4.54 X 10-^ 

0.01 

1 

2.15 


Table 1; GRB110426A parameter values. Values in bold text are not properly measured or reported and therefore 
default values are used ll52l . 



Figure 5: The fluences of the 57 selected GRBs (black 
curves and blue dashed curve) as generated by Neu- 
CosmA and their sum fluence (thick red curve). One 
GRB is brighter than the others by an order of mag¬ 
nitude above 10^ GeV (GRB 110426A, blue dashed 
curve). 


57 GRBs for the expected signal, and only searching in 
the 10 minute window surrounding each GRB. 

When optimizing our cuts, we use average, energy- 
dependent flavor ratios at Earth, which are calculated 
using the individual flavor ratios of each GRB in our 
sample, as output by NeuCosmA — the contribution 
of each GRB is weighted by its relative neutrino flu¬ 
ence. This is important, since electron neutrinos are 
more likely than other flavors to pass our trigger and 
analysis cuts due to charged-current events depositing 
the full neutrino energy in the particle shower. See Sec¬ 
tion |62] 

Among the set of analysis cuts described in the dif¬ 
fuse neutrino search, the Delay Difference Cut, the Re¬ 
construction Quality Cuts, and the Peak/Correlation Cut 
were re-optimized for this search. The three cuts that 
were re-optimized are all based on the quality of the 
directional reconstruction while the rest of the cuts are 
designed to reject specihc types of backgrounds such as 


CW and calibration pulser events. The Delay Difference 
Cut ensures that the reconstruction direction derived 
from all the borehole antennas of the same polarization 
is consistent with the delay observed between the sig¬ 
nals in the two antennas with the strongest signals. The 
Reconstruction Quality Cuts ensure that the event can 
be characterized by a single well-dehned pointing di¬ 
rection on the interferometric reconstruction map. The 
Peak/Correlation Cut requires that events have strong 
correlation between the signal strength and the cross¬ 
correlation value from the interferometric map, which 
is expected from impulsive events. 

A total of four cut parameters or options from these 
three cuts are allowed to vary to give the best expected 
limit on the dominant GRB event from the NeuCosmA 
model. For the Delay Difference Cut, we only consider 
whether to remove the cut, since it is largely redun¬ 
dant with other cuts. The Reconstruction Quality Cuts 
have two cut parameter values. Apeak Apeak/Atotab 
which ensure that the reconstruction direction is well- 
dehned and unique, respectively. Parameter Apeak is the 
maximum allowed area in square degrees on the inter¬ 
ferometric map surrounding the best reconstruction di¬ 
rection where the correlation remains high. Parame¬ 
ter Apeak/Atotai is the maximum allowed ratio between 
the high-correlation area around the best reconstruction 
direction and the high-correlation area from the entire 
map. The last parameter that was included in the opti¬ 
mization was the Peak/Correlation Cut Value, which is 
a unitless parameter that dehnes the minimum required 
value of a linear combination of the signal strength and 
the peak correlation value on the interferometric map. 

The expected number of neutrinos from each GRB 
and the background expectation based on the time of 
each GRB are obtained using the re-optimized cuts. For 
each GRB, we use its direction and predicted energy- 
dependent flavor ratio to obtain the analysis-level effec¬ 
tive area of the Testbed as a function of energy. The 
effective area of the f* GRB is obtained from 

the effective volume using the assumption that the di¬ 
mensions of the detector are signihcantly smaller than 
the interaction lengths lICTl : 


W) 


/int(F) 


( 2 ) 
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Cut 

Delay Difference Cut 

Reconstruction Quality Cut 

Peak/Correlation Cut 

Parameter 

On/Off 

-'^peak 

-^total-^peak 

Peak/Corr. Cut Value 

Diffuse Neutrino Search 

On 

< 50deg^ 

< 1.5 

> 8.8 

GRB Neutrino Search 

Off 

<140 deg"^ 

< 16.4 

> 7.6 


Table 2: Comparison of cut parameter values of the analysis. See text for details. 


where is the effective volume, calculated using 

Eq. Q, and /int(E) is the neutrino interaction length. 
The latter is given by 


lUE) 


mN 

^v~ice {E)p 

ice 


(3) 


where is the density of ice, o-y^icdE) is the 
cross-section of neutrino-nucleon interactions derived 
in Ref. ||68l , and m^v is the nucleon mass. 

The total expected number of neutrino events is 


^°p' = Z (/^ ’ 

(4) 

where i is the index of the GRB (total 57 GRBs) and 
F'(E) is the neutrino fluence [GeV 'cm“^] of the /* 
GRB . The factor In(lO) in Eq. Q is obtained by substi¬ 
tuting linear energy integration for logarithmic integra¬ 
tion, dEjE = £/ln(E) = In(lO) ■ c/log[Q(£’). 

Eigure shows the differential distribution of back¬ 
ground events as a function of the final Peak/Correlation 
cut. We estimate the expected number of background 
events by fitting an exponential function to this distribu¬ 
tion. 

As described at the beginning of the section, we de¬ 
rive the background estimate from the background anal¬ 
ysis window for each GRB, which is distinct from the 
signal window. We consider the background analysis 
window to be the hour on either side of each GRB time, 
minus the 10 minutes surrounding each GRB. The 55 
minutes on either side of a GRB (total 110 minutes) is a 
background analysis window and 5 minutes before and 
after the GRB is a neutrino signal window. A 10-minute 
period centered around the middle of the T’go window 
should be sufficient to encompass the expected emission 
period for all the GRBs examined in this study if we as¬ 
sume that gamma rays and neutrinos are produced si¬ 
multaneously. The 110-minute background period pro¬ 
vides sufficient statistics for a study of the background 
around the times of each GRB. This is the same method 
used in the ANITA GRB analysis ll34l . 

Using the data in the background analysis windows, 
we optimize our analysis cuts to give us the best ex¬ 
pected limit, and, using these optimized cuts, we obtain 



Peak/Correlation Cut Value 


Eigure 6; The differential distribution of events found 
in the background analysis windows of the 10% data 
set as a function of Peak/Correlation Cut Value after all 
other cuts have been applied. This distribution is fitted 
with an exponential function (red line) which is used to 
extrapolate the number of expected background events 
for a higher Peak/Correlation Cut Value. The optimized 
value is 7.6. 


the expected number of events from the background and 
signal windows. We compute the best expected 90% 
confidence level (C.L.) upper limit Eul on the neutrino 
fluence by minimizing 


Eul(E) = E,um(E) ■ ^ , (5) 

4 *exp 

where Esum(E) is the sum of the neutrino fluences from 
the 57 GRBs, Aexp is the expected number of neutrinos 
that pass the cuts, and Aul is the 90% C.L. upper limit 
on the number of signal events given the number of ex¬ 
pected background events. 

Table [^summarizes the final set of cut parameters af¬ 
ter the optimization. After the optimization, we expect 
0.072 events in the signal windows in the entire data 
set. This background expectation in the signal windows 
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is at approximately the same level as the expected back¬ 
ground events in the diffuse neutrino search, but now 
we achieve a factor of 2.4 improvement in the overall 
analysis cut efficiency for the summed fluence from the 
57 GRBs due to changing the analyzable time by a fac¬ 
tor of 566. To obtain the background expectations for 
the background windows in the 10% and 90% sets, we 
simply scale the 0.072 events by the livetime in each 
sample. In the background analysis windows in the 
10% subset, we expect 0.079 background events and no 
events survived. 

In the second stage of analysis, we look at the num¬ 
ber of events in the background analysis windows in the 
remaining 90% of the data set. This is to make sure that 
the background estimation derived from the 10% subset 
is consistent with what we see in the remaining 90% of 
the data. In the 57 GRB background analysis windows 
in the 90% data set we expected 0.72 events and two 
events survive. 

In the final stage of the analysis, we search in the en¬ 
tire data set for neutrino events in the signal windows 
surrounding the 57 GRBs over a total of 570 minutes. 
We used the same optimized analysis cuts defined in the 
first analysis stage. 

6. Results 

6.1. Upper limits on GRB neutrinos 

We expected 0.072 background events in the signal 
region in the entire data set and found no events. From 
NeuCosmA, the expected number of neutrino events 
from the 57 GRBs is 2.4 x 10“^. From simulation, the 
analysis efficiency for triggered events from the fluence 
calculated for GRB110426A is 6%. We placed a 90% 
C.L. limit on the combined fluence from the 57 GRBs. 

Figure [7] shows the total, or stacked, fluence from 
the 57 GRBs calculated with NeuCosmA, and the GRB 
neutrino fluence limit that we set from 10^ to 10^° GeV. 
At lower energies, the ARA Testbed sensitivity drops, 
and 10*° GeV is the maximum energy with which Neu¬ 
CosmA emits neutrinos. 

In order to compare our limit with those from other 
experiments that used a different set of GRBs for their 
analyses, we also provide the inferred quasi-diffuse all¬ 
flavor neutrino flux limit. This assumes that the average 
fluence of the 57 analyzed GRBs is representative of 
the average fluence from GRBs for any other extended 
period. With this assumption, the quasi-diffuse neutrino 
flux limit is 

1 V° 

F^d) = E^F X — , (6) 

An Vgrb 


10 ® 



Neutrino Energy [log^g(E/GeV)] 


Figure 7; The 90% upper limit on the UHE GRB all¬ 
flavor neutrino fluence from 57 GRBs. Total fluence 
from NeuCosmA for the 57 GRBs is shown with a red 
shaded area and the limit from the ARA Testbed above 
10^ GeV is shown with a black solid curve. 

where E^F is the fluence limit, Vgrb = 57 is the num¬ 
ber of analyzed GRBs, and is the average num¬ 
ber of GRBs that are potentially observable by satellites 
per unit time llTSll . and is chosen as 667/year to be con¬ 
sistent with the IceCube and ANTARES GRB neutrino 
searches Il3^l3^ . 

Eigure shows the quasi-diffuse neutrino flux limit 
from ARA and other experiments. Our limit is the first 
UHE GRB neutrino quasi-diffuse flux limit at energies 
above 10’ GeV. The sensitivity of IceCube extends to 
this energy region, but their quasi-diffuse limit is pub¬ 
lished only below 10’ GeV, where their sensitivity is 
greatest. 

6.2. Ejfects of uncertainties and model parameters 

Our calculations are unavoidably affected by uncer¬ 
tainties in the values of astrophysical parameters — on 
which we expand below — and of particle-physics pa¬ 
rameters, including cross sections, multiplicities, and 
lepton mixing parameters. Astrophysical uncertainties 
affect each source in a different way, and, in a source 
sample, partially average out. Particle-physics uncer¬ 
tainties systematically affect the fluxes from all sources 
in the same way, but are considerably smaller than as¬ 
trophysical uncertainties; see, e.g., Eig. 19 in Ref. m 
for the effect of the uncertainty on the mixing parame¬ 
ters. We have therefore assumed in our calculations the 
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central values of the particle-physics parameters. 

In the calculation of our limits, we assumed nominal 
values of the astrophysical model parameters. We now 
comment on the effect of varying these values. Ref. 1721 
showed the effect on the shape and flavor composition 
of the diffuse GRB neutrino flux of assuming distri¬ 
butions of values for the magnetic field intensity, bulk 
Lorentz factor, and shape of the source photon spec¬ 
trum. In stacking analyses, the combined uncertainties 
on astrophysical model parameters can lower or raise 
the quasi-diffuse flux by one order of magnitude Il43]l . 
The baryonic loading is particularly poorly known; in 
our analysis, we adopted the commonly used value of 
10 for all bursts l43l . In reality, it could be lower 
or higher by a factor of 10. Since the baryonic loading 
linearly scales the neutrino flux, this would shift the flux 
down or up by one order of magnitude ll4^ . 

Another source of uncertainty is the finite size of the 
GRB sample used to derive the quasi-diffuse flux. For 
instance, the uncertainty associated to the discrete sam¬ 
pling of the underlying redshift distribution of GRBs 
ranges from 56%-72%, for a sample of 50 bursts (the 
present analysis uses 57 bursts), to 25%-28%, for a 
sample of 1000 bursts (90% C.L.) Il72l . 

While we have considered GRB jets whose baryonic 
content is dominated by protons, GRBs might be able 
to synthesize 1173117^1751 and accelerate 117^177117811791 
l59l [60l nuclei. If nuclei can reach energies as high as 
protons, neutrino fluxes are comparable ll59l : otherwise, 
neutrino yields from nuclei could be up to two orders of 
magnitude lower llTTl . An exploration of GRB neutrino 
limits assuming different jet mass compositions is be¬ 
yond the scope of this paper. 

Alternative fireball emission models, such as sub- 
photospheric lISOl ISTl l82l l83l l84l l85l and magnetic re¬ 
connection ||8^1871|8811M1 models, may result in quasi- 
diffuse neutrino fluxes up to one order of magnitude 
lower than the flux from the internal-collision model we 
adopted Il371l38l . 

While our results in Figs.|7]and|^use average, energy- 
dependent flavor ratios at Earth (see Section 5.2 1 , we 
considered the impact of variations in flavor ratios. In 
Ref. ||90l , it is argued that for (1 : 2 : 0)s flavor ratios 
at the source, high-energy neutrinos from astrophysical 
sources can reach Earth with ratios (r ; 1 ; Ije where 
0.57 < X < 2.5, and Ref. HI] finds an electron frac¬ 
tion between 20% and 59%, corresponding to the range 
0.5 < X < 2.9. Eor (1 ; 1 ; Ije ratios in the incident 
flux, at the trigger level the ratios of detected neutrinos 
become (2 ; 1 ; 1), and, at the analysis level, they be¬ 
come (6 ; 1 ; 1). Due to this effect, neutrino fluxes 
with flavor ratios of (0.5 : 1 : Ije and (2.9 : 1 : Ije, 


with the same all-flavor normalization, would result in a 
25% lower and 50% higher number of neutrinos passing 
the trigger and analysis cuts, respectively, and a corre¬ 
sponding weakening or strengthening of the limits. 

7. Future prospects 

Eor future analyses using two ARA deep stations, we 
expect to have at least a factor of 6 improvement in 
sensitivity compared to this one using Testbed data as¬ 
suming the same analysis with similar cuts. There is 
a factor of ~ 3 expected increase going from the shal¬ 
low Testbed station to a 200 m deep-station and another 
factor of ~ 2 for the number of deep stations currently 
operating. In addition, we plan to increase the num¬ 
ber of deep stations. Eig.[^ shows the expected ARA37 
trigger-level limit based on these and other improve¬ 
ment factors similar to those described for the diffuse 
neutrino search ll29l . Below, we motivate an expecta¬ 
tion for a high analysis efficiency in future ARA GRB 
analyses. Eurthermore, the implementation of a phased 
array trigger design, as described in Ref. Il92l . currently 
funded for an initial deployment in 2017-2018, would 
decrease the trigger threshold and improve the sensitiv¬ 
ity to neutrinos from GRBs. 

In the future, by restricting our GRB searches in di¬ 
rection (so as not to include the South Pole direction), 
and by improving the way we reject CW backgrounds, 
we expect that we may eliminate all cuts but those de¬ 
signed to reject thermal noise. ARA has the ability to 
reconstruct the directions of RE signals, and we plan to 
develop the capability of reconstructing neutrino direc¬ 
tions also, using polarization and spectral information. 
In addition, we are working to replace our CW cuts with 
biters. Keeping only cuts designed to reject thermal 
noise would leave the Reconstruction Quality Cut and 
the Peak/Correlation Cut as those with an important im¬ 
pact on our sensitivity. With only these cuts, we bnd 
that the analysis efficiency for the dominant GRB bu- 
ence in this paper increases from 6% to 14%, a factor of 
2.3 increase beyond the increases mentioned above due 
to expansion of the array. 

Improvements in the reconstruction by using an algo¬ 
rithm that solves for event distance and additional anten¬ 
nas in design stations are expected to lead to improve¬ 
ments in the analysis efficiency by an additional fac¬ 
tor of a few. Although the Reconsttuction Quality Cut 
was signibcantly relaxed here compared to the diffuse 
analysis ll2^ . its efficiency against simulated triggered 
events was ~ 30%, primarily rejecting events with a low 
signal-to-noise ratio (SNR). Improvements to the recon¬ 
struction method under development will be able to in- 
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crease the efficiency of reconstructing these low-SNR 
events. Additionally, in the design stations, the number 
of pairs of antennas of each polarization contributing to 
the interferometric map increases from 6 to 28, which 
is expected to improve the efficiency, in particular, by 
giving low-SNR events a higher peak correlation value 
to differentiate it from noise. 

8. Conclusions 

Using data from the ARA Testbed station from Jan¬ 
uary 2011 to December 2012, we have searched for 
UHE neutrinos from GRBs. We selected 57 GRBs that 
occurred during this period within the field of view of 
the Testbed. We searched for GRB neutrinos in a time 
window around each burst. The resulting reduced back¬ 
ground allowed us to loosen our analysis cuts and im¬ 
prove our analysis efficiency for neutrinos from the 57 
GRBs by a factor of 2.4. The GRB neutrino spectra 
were calculated using NeuCosmA, an advanced high- 
energy astrophysical neutrino fluence generator. 

We found zero events passing the cuts for our search, 
which is consistent with the expectation. We obtained 
a GRB neutrino fluence limit and the first quasi-diffuse 
GRB neutrino flux limit for energies above 10^ GeV. 

Future analyses from two ARA deep stations are ex¬ 
pected to have at least a factor-of-6 improvement in sen¬ 
sitivity compared to the present analysis with the ARA 
Testbed, assuming the same cuts. Another factor of 
about 10 is feasible from planned developments in re¬ 
construction and CW filtering capabilities at the analy¬ 
sis stage with the current deep station design. 
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Figure 8: The ARA-Testbed quasi-diffuse all-flavor flux limit. We include limits from IceCube lIMI and 
ANTARES ll^ for comparison; we have multiplied them by a factor of 3 to make them all-flavor. IceCube re¬ 
cently published a search for neutrinos from GRBs based on four years of data ED, but did not include a limit on the 
quasi-diffuse flux. Preliminary estimates indicate that the latest result would improve upon the IC40-I-59 limit shown 
here by about an order of magnitude. The ARA37 limit is the trigger-level sensitivity based on scaling the Testbed 
using factors described in the diffuse neutrino search Il29l . For reference, several diffuse limits have been included (in 
grey); the Testbed diffuse flux limit Il29l . the ARA 2-station diffuse limit 1301 . and the 2012 Extremely High Energy 
(EHE) diffuse limits from IceCube 1691 . The points in grey represent the fluxes from the IceCube high-energy starting 
events (HESE) using 3 years of IceCube data CQl. For comparison, the Waxman-Bahcall upper bound on the neutrino 
flux from UHECR thin sources is 3.4 x 10“® GeV cm“ 
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